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Abstract. We diagonalize the XX model with a finite number of spins and periodic 
' (— I ', boundary conditions. We solve for the ground state, focus on the rapidity of the 

convergence to the thermodynamic limit and study the features of multipartite 
entanglement. 

^ ' 1 Introduction 

Recent advances in quantum information processing motivated a widespread interest for spin 
systems as models for quantum computers [1]. Spin chains are often studied, in particular in 
^ , the thermodynamic limit [2|3|4j . where one can neglect the contributions deriving from the 

■ finite size of the system, that scale like 0{1/N), with N the number of spins. However, due 
Q>^ I to experimental and theoretical difficulties, such as decoherence and/or imperfections in the 

' quantum hardware, nowadays it is only possible to assemble and control the interaction of a 

small number of qubits (always less than 10). For this reason, it is extremely important to take 
into account all finite-size effects in the system Hamiltonian. 
• The investigation of the last few years has focused on entanglement [5] in diverse finite- 

OO ! size models, by means of direct diagonalization [6]. These studies were boosted by the recent 

■ discovery that entanglement can detect the presence of quantum phase transitions [7]. In this 
article we will study the finite-size ID quantum XX model with periodic boundary conditions. 
We will exactly diagonalize the Hamiltonian by a Jordan- Wigner transformation followed by a 
deformed Fourier transform. Some of the results wc will obtain were alluded to in the seminal 

JH I article by Licb et al. ^ and were discussed by Stclmachovic and Buzek [Hj and very recently by 

Canosa and Rossignoli [9] . We will focus here on the features of the ground state and evaluate 
the "forerunners" of the critical points, as the quantum phase transition (QPT) occurs only 
in the thermodynamical limit. By introducing a finite-size parameter, we will show that the 
thermodynamical limit is very rapidly approached, so that a 10 qubit system already represents 
a very good approximation. Finally, we shall analyze the multipartite entanglement features of 
the finite-size ground state. 

2 Exact diagonalization of the finite size XX model 
2.1 Exact diagonalization of the Hamiltonian 

The XX model for a collection of N qubits on a chain is described by the Hamiltonian 



i=0 




Fig. 1. A spin ciiain witli periodic boundary conditions. 

where J is a constant with dimensions of energy and g a dimensionlcss parameter proportional 
to the transverse magnetic field. We consider periodic boundary conditions (cyclic chain, see 
Fig. [ID 

CTo — (TN. (2) 




The usual diagonalization scheme is a two-step procedure based on the Jordan- Wigner ( JW) 
transformation followed by a Fourier transform. The use of the JW transformation is based on 
the observation that there exists a unitary mapping between the Hilbert space of TV qubits 
(C^)®^ and the antisymmetric Fock space JF_(C^) of spinless fcrmions on N orbitals. We can 
define annihilation and creation operators for JF_, 

Vi e Zjv, (3) 

where = {af iticrf )/2, Zyy = {0, 1, . . . , A'^ — 1}, and is the number operator counting the 
holes between and i — 1, 

"a-E(i--K)=E-7-;- (4) 

j=0 j=0 

It can be shown that the JW operators satisfy anticommutation relations 

{c„cj} = 0, {clc]} = 0, {c„c]} = %, Vi,je {0,1,...,7V-1}, (5) 

whereas the corresponding Pauli operators anticommute only on site: 

K±,(T±} = for i = j, [at,af] = for t ^ j. (6) 

The periodic boundary conditions assigned to the cyclic chain do not hold in general for the 
fermionic operators indeed, by prolonging ([31) up to z = TV we would obtain 

co = fT(7, cjv = e'™ia^ = e'™i(7o- =e'™ico, (7) 

where = iini is the number operator that counts the total number of spins down (and the 
number of holes), from to TV — 1. As a consequence, the Hamiltonian written in terms of the 
JW operators is characterized by the presence of a boundary term: 



H = -.] 



N~l N~2 



J2 5(1 - 2c,c]) + J2 + + (e"("^+'^ - l)(cjv-i4 + coc]v„i) 



(8) 
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In the thermodynamic limit [3"4] the XX Hamiltonian is diagonalized by introducing the discrete 
Fourier transform 

— ^e^Cfc, VjGZat (9) 



and neglecting the last term in ([8]), since its contribution scales like 0{1/N). We intend to 
evaluate the corrections due to the finite size of the chain. 

The main difficulty introduced by the boundary term in the Hamiltonian ([8]) is that it breaks 
the periodicity of the JW operators, due to the arbitrary dependence of the phase e"^"'^ on the 
ordering of the spin on the circle. This phase clearly depends on the state the Hamiltonian H 
is applied to. However, notice that the parity of the spin-down number operator is conserved 
(although not so the number operator itself) 

[e'™Si7]=0. (10) 
Therefore, let us consider the parity operator 

V = e'''("i+i) (11) 
and the spectral decomposition of the number operator n^, 

N 

ni = 51 '^il"i>("il' (12) 

114^=0 

where |n|) is the eigenstate of nj^ belonging to the eigenvalue ni . In this basis the parity operator 
can be written as: 

N 

where 

p+= P-= KX^il (14) 

odd ni even 

are the projection operators associated to the eigenvalues ±1 of respectively. Due to Eq. 
(|10p . the Hamiltonian preserves the parity sectors and can be decomposed as 

H = P+HP+ + P^HP^ = + H^'-l (15) 

The analysis can then be separately performed in the two parity sectors, where V acts as a 
c-number. 

Let us define a deformed Fourier transform 

N-l 

cj = 7=e " « Cfc, (16) 

where we added a position dependent phase aj (j G Zjv denoting the site). The anticommutation 
relations for Ck still hold: 

{cfc,4'} = 0, {4„4,} = 0, {ck,cl} = Skk'. (17) 

The phase e » can be determined by imposing that the last term of Eq. ([8|), after Fourier 
transform, has the same form of the other A'^ — 1 terms: namely we require that the following 
two expressions, 
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and 



'W-l _27riQo 1 \ ^ \ ^ 2xi(W-l)fc 



Ckcl 



fe=0 fc'=0 



gi7r(rn + l)g 



JV-l _27riao i \ ^ \ ^ 2iri(W-l)fc 2wiNk' . „f 

" ]^ !L 2^ ^ " ^ " ^^^fc' 

fc=0 fc'=0 



have the same phase. It follows that 



e iv = e'^'^^i+i^e 5v . 



The solution is given by 
where a satisfies the equation 



a, = ja + ao, 



(19) 

(20) 
(21) 
(22) 



On the other hand, the phase of the first site ao is a free parameter. The two solutions in the 
two parity sectors are 



_ I mod N if V = +1 {n^ odd), 
I i mod if V = —1 {rii even). 

Summarizing, by using the (sector dependent) deformed Fourier transform 

JV-l 

Cj = —-=e jv e « > e « Cfc 
'^-^ fc=0 

and substituting into Eq. ([8| we obtain 

N~l , 

H=-2JY, etc,-- 

fc=0 ^ 



(23) 



(24) 



g — cos 27r 



(25) 



where a is given by (|23p . In this way, if is diagonalized and the boundary term removed. 
Therefore, the Hamiltonian is given by (|15p with 



N-l 



fe=0 ^ 
N-l , 



g — cos 27r 



k=0 



q — cos 27r 1 



(26) 
(27) 



2.2 Energy spectrum 

In this section we will focus on the spectrum of the system. We shall set henceforth J = 1. The 
state with no fermions has an energy density that is independent of a and of the number of 
sites N. Indeed, from Eq. (pS)) . 



N N ^ 

k=0 



COS 2tt 



a + k 
N 



(28) 
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If we add one fcrmion the energy density reads 



£i{k,g) = evacCg) 



where a depends on N: 



2 



g — cos 27r 



a 



N 



(29) 



N even 
A^odd : 



iV - 1 odd 
— 1 even 



a = mod 
a = i mod N. 



(30) 



In Fig. [2] we represent the single particle energy spectra corresponding to iV = 8 and 9 sites 
(representative of an even/odd number of fermionic sites, respectively). Different curves are 
parametrized by A: G Z^v and one notices the presence of degeneracies in both cases. Since 
we are interested in the ground state of the system, we focus on the lowest energy levels and 
consider the values assumed by the function cos[27r(Q; + k)/N] in the four possible cases [2 
parities of number of sites N and 2 values of a as in Eq. as shown in Fig. [31 Notice that 

these results can be described in terms of regular polygons inscribed into a circle of unit radius, 
see Fig. |4l 



N even (iV = 8) N odd {N = 9) 




Fig. 2. Single particle spectra ei{k,g) (solid lines) and vacuum density energies evac{g) (dashed lines). 
Different hues correspond to different k £ Zn according to 1)29^ . 





Fig. 3. Plot of cos[27r(Q + k)/N], k € Zn for A/' = 8, 9 and q = mod N, i mod N. 
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From Fig. [3] one obtains the values of k that minimize the energy per site; in the 1-particle 
sector one has 

f TV even =^ a = mod TV ^ fc = ^, 

\ (31) 
[A^odd => a = i mod TV ^ fc = 

Similarly, in the 2-particle sector the minimum energy is at 

=> a = i mod TV ^ {/ci,fc2} = {f - 1, f } , 

a = mod N {ki,k2}= {^,^} . ^^^^ 

As a result, the general expression of the lowest energy levels in the different n-particlc sectors 
does not depend on the parity of TV. For n fermions, one gets 

,33, 

In Fig. O we plot the lowest energy levels corresponding to < n < TV for TV = 8 sites. The 
intersections of levels corresponding to n and n + I fermions (starting from n = 0) define 
the level crossing points or quantum critical points Qc, where an excited level and the ground 
state are interchanged. The analytic expression of the critical points is easily obtained by the 
condition er"(gc) = e"|"i(ffc). We find 

sin(WTV)-sirr[(n+lWTV] ^^^^^^ 
' sin(7r/TV) , _ _ \ ) 

As a consequence, the ground-state energy per site is 

/ X / 2n\ 2 sin(7i7r/TV) . , , , ^, , ,. . 

^^h^ sinWV 5e(.g.(n-l),5c(n)), (35) 

with < ?i < TV + 1 and .gc(-l) = -oo, gc{N + 1) = +oo. Thus, for g e [gdn - l),gc{n)), the 
ground state contains n fermions. Note that gdO) = —1 and gc{N) = +1, independently on TV. 

In the thermodynamic limit the critical points correspond to the QPT. In order to study 
the properties of the ground state in this limit, we evaluate the envelope of the lowest-energy 
levels. This is obtained by the equation 

5e™'"fo) TV 

" - ^ n(.g) = -arccos(-5X^), (36) 
on TT 

where 

_ sin(^/A^) 

TT/TV 
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A^ = 8 




Fig. 5. Lowest-energy levels e™'"((;) for different number of fermions n; the intersection between the 
energy levels corresponding to n and n + 1 fermions (starting from n = 0) are the quantum critical 
points (dots). 



is a finite-size parameter depending on the number N of sites in the chain. The envelope is 
obtained by plugging n{g) from Eq. ([36]) into Eq. (|35p . By noting that n = for g < — I/xat 
while n = N for g > 1/xn, one gets 



2 , A 2 VT 

g [1 arccos(-gxw) 



for I5I < — , 

XN 



for I5I > 



XN 



It is easy to see that the envelope and its first derivative arc continuous functions at \g\ 
whereas the second derivative 



XN 



d'ee 



for I5I < — , 

XN 



dg2 



(38) 



(39) 







for \g\ > 



XN 



diverges to —00, for \g\ f 1/xn- 

In Fig. [HI we plot the lowest-energy levels and the ground-state envelopes for iV = 9 and 45 
sites. Notice that the scales of the energy and magnetic field are similar. This result is confirmed 
in Fig. [71 where we plot the envelopes for = 9 and 45; their difference is negligible if compared 
to the envelope for 5, which shows that = 9 is already a good approximation of the 

thermodynamical limit. 

It is interesting to evaluate the finite-size effects. This can be done by considering the 
dependence of xn on N. We plot this function in Fig. [5] (left). For N > 10, xn — 1 with 
excellent approximation. We can obtain further information by looking at the thermodynamic 
limit of the envelope ((38|l , 



lim ( 

Af— >oo 



1 arccos(— g) 



2 r- 

n 



\9\ 



g2 for \g\ < 1, 

for |g|>l, 



(40) 



which is nothing but the ground state energy density of the thermodynamic system. The relative 
error reads (for g = 0) 



Ae 



.(O)-goo(O) 

koo(o)| 



1 

Xn 



1 



1 

3! \N 



for — > 00. 



(41) 
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AT = 9 iV = 45 




Fig. 6. Ground-state envelopes £cnv{g) for = 9 and 45 sites. 
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Fig. 7. Ground-state envelopes £cnv(i7) for A'' = 5, 9, and 45 sites. 

iV 



1.0 1 ^ 0.00 




10 20 30 40 50 



TV 

Fig. 8. Finite-size parameter Xiv (left) and relative error zi£/|e| (right) for A'^ = 1, . . . , 50. 



It is interesting to notice that the relative error is OiXjN"^), which is better than what one 
naively expects, see comment after Eq. ([9]). 

Figure[8] (right) confirms this result. The effect of the finite size of the chain can be neglected 
already for a relatively small number of sites. 
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2.3 Ground state revisited 



It is interesting to reinterpret the ground state of the chain in terms of spins. For g < —1 the 
ground state is the state with no fermions. By definition, this corresponds to all spins down 
(antiparallel to the z axis) in the chain 



\i'o) = m)---\i)- 



(42) 



This means that for strong magnetic fields the interaction between the field and the spins is 
larger than that among the spins in the xy plane. 

As the magnetic field tends to zero, the number of spins in the down state decreases and 
the general expression of the ground state for g G [gc{n — ^)igc{n-)] is given by 



1 

N 



E 



i\<n<-<3n 



U)oU)i---U)iv-i 



N 



E [^Jl J2,...J,.(-1)"^'(-1)'""^^'^'"^'H-1)'"~^^'^'"^'-' ■ • ■ (-1)^""^"-^ 



J1<J2<---<J„ 



X U)o---|T),,---|T),2---ITk---U)^-i, 

(43) 



where ^ji,j2,---,jn the sum over all permutations {1, 2, . . . , n} — s- {pi,p2, . . . ,pn} 



E 



(_]^)Pg^ C^iJpi +k2jp2 H hk„jp^ ) 



(44) 



Note that in this case the state of the system cannot be written as a tensor product and {tpn) 
is entangled. When the magnetic field is very strong and positive g > 1 (parallel to the z axis) 
all spins are aligned along its direction and the state of the chain is 



l^^> = IT)IT>---|T)- (45) 

These properties are summarized in Fig. O We will analyze the entanglement structure for a 
finite-size chain in the following section. 
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Fig. 9. Classification of the ground states in terms of entanglement {N = 8). 
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3 Multipartite entanglement of the ground state 

3.1 Probability density function characterization of multipartite entanglement 

The evaluation of the entanglement stored in the ground states of quantum spin chains has 
become a central problem in this field of research. The study of the entanglement properties of 
the ground state usually explores the link between quantum phase transitions and entanglement 
generation and previous research mainly focused on the thermodynamic limit [7]. 

One can use several measures of entanglement [10] . A typical approach consists in considering 
a given bipartition and evaluating the bipartite entanglement (in terms of the von Neumann 
entropy, the linear entropy, the purity, or in general a well defined entanglement monotone 
[5]). A possible approach to multipartite entanglement is to analyze the statistical properties 
of bipartite entanglement over all balanced bipartitions (each part consisting of one half of 
the chain) [llj . This characterization has been applied to the ID quantum Ising model in a 
transverse field, yielding interesting results |12j . 

We will characterize in a similar way the multipartite entanglement of the XX chain. Con- 
sider a chain of N spins and consider a partition in two subsystems A and B, made up of 
Na and Nb qubits {Na + Nb = N) , respectively. For definiteness we assume Na < Nb ■ The 
total Hilbert space is the tensor product Ti, = Ha ^ "Hb with dimensions dimTi^ = 2^-*, 
dimHB = 2^«, and diuiH = 2^^+^« = 2^. 

Let us denote the ground state by lijjgs) and consider the purity of subsystem A (which 
equals that of subsystem B) 

Ti-AsdV'gs)) Tr^p^, PA = Trap, p = |?/'gs)(V'gs|, (46) 

Tta (Trs) being the partial trace over subsystem A (B), and take it as a measure of the 
bipartite entanglement between A and B. We note that 

2-^^ < TTAB < 1, (47) 

where the minimum (maximum) value is obtained for a completely mixed (pure) state pA- 
Therefore, a smaller value of ttab corresponds to a more entangled bipartition {A,B). ttab 
will depend on the bipartition and entanglement will be distributed in a different way among 
all possible bipartitions. The key idea is that the average of t:ab, denoted by p, measures the 
amount of entanglement (the smaller p, the larger the entanglement), whereas the standard 
deviation a measures how well this entanglement is distributed among bipartitions. Clearly, if 
the distribution function ttab is not very regular, higher moments will be necessary in order to 
properly characterize it. 

3.2 Average and standard deviation of the distribution 

We numerically evaluated the distribution of bipartite entanglement for the finite-size XX model 
on the circle. In Fig. [TO] we plot the average entanglement p over balanced bipartitions versus 
the coupling g (for iV = 4, . . . , 10 sites). 

Note that, as expected, the maximum entanglement is reached at g = 0. Moreover, the 
larger the system, the smaller p. An interesting peculiarity of the XX model, at variance with 
the Ising chain [12] . is that p{g) is not continuous: we observe jumps between plateaux. These 
jumps correspond to the level crossing points defined in Eq. (|34p . In particular, for \g\ > 1, the 
ground state of the model is completely factorized and, consequently, we have no entanglement 
{p ~ 1). The largest jump is found at |g| = 1, where p ~ 0.5 for all TV. 

The standard deviation of the distribution is plotted in Fig. [TT] We displayed separately the 
cases N even (left) and odd (right). Also in this case we observe a stepwise behavior, depending 
on the number of fermions in the ground state. For odd A'^ the maximum of a tends to decrease 
for a larger chain: the curves are monotone for — 1 < g < and < g < 1. On the other 
hand, for even N, a shows a more complicated structure; the curves are not monotone. The 
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Fig. 10. (Color online) Average fi of ttab over all balanced bipartitions for the XX chain with periodic 
boundary conditions, with A'' = 4, . . . , 10. 



a 

O.lOr 



0.08 - 

ofofT- 



0.04 
0.02 

-1.0 -0.5 0.0 



iV = 4 



iV = 6 



Af = 8 
JV = 10 



^,9 



a 



0.10 
0.08 



0.04 
0.02 

-1.0 -0.5 0.0 



Af = 5 
N = 7 
N = Q 



:,9 



Fig. 11. (Color online) Standard deviation a of ttab over all balanced bipartitions for the XX chain 
with periodic boundary conditions, with an even (left) and an odd (right) number of sites. (Left) 
N = i,6, 8, 10. (Right) 'n = 5, 7, 9. 



maximum is not a decreasing funetion of the size (for = 8 it is larger than for = 6). As 
a general trend, however, for both even and odd N, a tends to decrease with A^. Additional 
investigation is required in order to understand the features of multipartite entanglement and 
in particular the large- A^ limit. As explained in [T^], different scenarios are possible, depending 
on the behavior of a/ii: if this ratio tends to zero with A^, entanglement will tend to be well 
distributed among different bipartitions. 



4 Conclusions 



We analyzed the finite-size XX model on the circle. We diagonalizcd the Hamiltonian by defining 
a deformed Fourier transform that takes into account the periodic boundary conditions. We 
also studied the ground state of the chain, by focusing on the points of level crossings, that 
forerun the QPT in the thermodynamic limit, and defined a finite-size parameter that quantifies 
how rapidly this limit is approached. Finally, we looked at the properties of the multipartite 
entanglement of the ground state in terms of the distribution of bipartite entanglement. There 
is considerable interest in the study of entanglement for quantum spin chains, both in view of 
applications and because of their fundamental interest. Future activity will focus on the study 
of different models, with more general interactions and/or boundary conditions. 
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